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Q\ i ABSTRACT 



We discuss the role of nonfactorized contributions in B — > ip + K* and charmed 



, meson decays. We demonstrate, using Dg — > tfiir as a model, how nonfactor- 



ization, annihilation and final state interactions can be built into effective and 
unitarized a± and 02- 

1. Introduction 

Factorization approximation, where the matrix element of a product of two color- 
singlet weak currents is approximated by the product of the matrix elements of 
£Cq 1 the individual currents, is the most commonly used scheme to calculate the am- 

plitudes for two-body hadronic decays of B and D mesons @. It was shown recently i 
that this approximation failed to reproduce the longitudinal polarization observed in 



B — > ip + K* decay in all the commonly used models of form factors . Subsequently 
Qh| it was shown i how a small nonfactorization contribution helps in understanding the 

longitudinal polarization in B — > ip + K* decay in all the commonly used models of 
form factors . 

In the following I will discuss three topics : (I) Nonfactorizatation and polarization 
in color-suppressed B — > ip + K* decay, (II) Nonfactorization in Cabibbo-favored D 
decays and (III) Nature of aiand a 2 . 



2. B -»• if) + K* decay 

The effective weak Hamiltonian for decays of kind b — > sec is : 

H w e ff = ^V cb V c * s {d(c6)(sc) + C 2 (cc)(sb)} (1) 

where (cb) etc. represents color-singlet (V-A) current, and V cs are the relevant 
Cabibbo-Kobayashi-Maskawa (CKM) mixing parameters and C\ and C2 are the QCD 
coefficients for which we use the values i: 

d = 1.12 ±0.01, C 2 = -0.27 ±0.03 (2) 



If we Fierz-transform in color space, we can write 

(cc)(sb) = ±(cb)(sc) + W (cX a b)(sX a c) (3a) 

and 

(cb)(sc) = ^(cc)(sb) + W (c\ a c)(s\%) (36) 

where A a are the Gell-Mann matrices and N c = 3 is the number of colors. It is conve- 
nient to define two parameters a\ and a 2 as follows: 

ai =d + -J-C 2 = 1.03 ± 0.01 (4a) 

a 2 = C 2 + -J-Ci = 0.10 ±0.03 (46) 
The decay amplitude for B —>■ ip + if* can be written in the following form : 

A(B -> ^ + if*) = -|K*V£ {a 2 < if*V> I (cc)(s6) I 5 > +Ci < K*if) \ H$ \ B >} 

(5) 

where 

^ = ^E(^ a c)( S -A a 6) (6) 

is a product of color-octet currents. In the factorization approximation, only the first 
term in Eq.(5) is kept. 

Confining ourselves to the factorization approximation, we can write the following 
expressions for the decay amplitude in terms of the relevant form factors : 

A(B^ij + K*)= Gp^m^f^ I (m B + m K ,) A^* (m 2 ,) 

2 

;ei ■ Vb£-2 ■ PbA 2 BK * (m^ 

n K * i 

2i 



(m B + m K *) v 

-^ pa ex^2 V VK* P VB°V BK \m%) \ (7) 



(m B + m K * 

where Gp — ^V c bV* s and t\ and e 2 are the polarization vector of the ip and K* 
respectively.From Eq.(7) one obtains! 

B{B^i, + K*) = 2.84a 2 2 | A x BK \m%) | 2 (S LL + S TT )% (8) 

where for longitudinal and transverse states of polarization, 

S LL = (a - bx)\ S TT = 2(1 + cV) (9) 



Figure 1: The domain of x and y allowed by the polarization data Eq.(12). Points 
A,B, . . . represent predicted values of x and y in various models. See Ref. [2] for 
details. 



x = tS ' y = \ ii5k ( 10a ) 



with 

M{m%) = V(m 2 ; 

and 

a = 3.147, b = 1.297, c = 0.434 (106) 
From Eq.(9), we obtain the longitudinal polarization as follows: 

\2 



bxY 



(11) 



S LL + Yjtt (a - bxf + 2(1 + cV) 
Experimentally S@ 

P L = 0.78 ±0.07. (12) 

In Figure 1 we have shown the allowed domain in (x-y) plane to one standard 
deviation of the central value of Pl data. We have also shown the value of x and y 
predicated by six different models. For details the reader is referred to Ref. [2]. 

It is clear that the commonly used models of form factors generate points in 
(x-y) plane that are several standard deviations removed from those required by 
polarization data. 

Assume now that there are nonfactorized contributions to the decay amplitude 
of Eq. (5). These could arise from two sources : (i) nonfactorized contribution to the 

first term comprised of color-singlet currents and (ii) contribution from H w which is 
a product of two color-octet currents. This latter contribution is enhanced relative to 
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Figure 2: Domains of x and y allowed by polarization data for different values of the 
nonfactorization parameter x- 



the former by the fact that | C\ |~ 10 | a 2 |. If for simplicity we assign nonfactorized 
contribution to the Lorentz structure belonging to A\ only, then we need to modify 
our formalism with the replacement 

A ± BK * (ml) -> A 1 BK * (ml) + + — Ai (8)n/ (13) 

where A^ 11 ? and A^ 11 -^ are the two nonfactorized contributions referred to above. 
The resulting formula for Pl becomes 

Pl = ( 4 ~ hX? (14) 

K - bx) 2 + 2 (e + cV) 

where 

f = 1 + — X (15a) 

X = (^i (8)n/ + ^i (1)nf ) /Ai BK *(ml) (156) 

In Figure 2 we have plotted P L for different values of x which is a measure of 
nonfactorized contribution . 

We note that for x ~ 0.1 all models of form factors accomodate data. We remark 
that if nonfactorization were put in, say A 2 , we would have needed much larger values 
of the parameter analogous to x- Thus even as little as a 10% nonfactorized contribu- 
tion provides an understanding of the longitudinal polarization in all commonly used 



and 



models of form factors 



3. Nonfactorization in Cabibbo-favored D decays. 

In the following we will discuss the role of nonfactorization in a selected few 
Cabibbo-favored D decays: D — > ttK,ttK*, pK. 
The relevant effective Hamiltonian is 

H w eff = ^VcsV u * d {C^ud^sc) + C 2 {uc)(sd)} (16) 

where for C\ and C 2 at charm mass scale we adopt the following values H@ , 

C l = 1.26 ±0.04, C 2 = -0.51 ± 0.05 (17) 

so that with N c = 3, we obtain 

at = 1.09 ±0.04, a 2 = -0.09 ±0.05 (18) 

Fierz transformation of the current products in color-space leads to 

(uc)(sd) = ^-(ud)(sc) ± ^ E (u\ a d)(s\ a c) (19a) 

and 

(«d) (sc) = — (uc) (sd) + ^Y1 (« A<lc ) (^ arf ) ( 196 ) 

H w {8) = \j2 (u\ a d) (^ a c) ( 20a ) 

* a 

Using the above equations we calculate the decay amplitude for D — > Kit, K*ti and 
Kp decays. As an example, 

A(D° -> IT**) = ^Vc^Cmi, - m^)a 1 e ^F ^K) (21a) 
A(D° - ^V) = ^K S KV*K> - mlW^F^(m\) (216) 



Define 



and 



and 



where 



1 F DX (m2) ai F DK (ml)j 



and 

a 2 «// = a 2 1 + „° 9x + ^ „°, 9x }> (226) 

Fq^ 71 -? and Fq 1 ^^ are the nonfactorized contributions from the product of color- 
singlet currents and F ^ n ^ and F ^ * arise from and Hw^ of Eq. (20), re- 
spectively. 

^From Eq. (21) we calculate the isospin amplitudes Bi, Ax and As, by setting 
the phases equal to zero in the following isospin decomposition, 

A(D° -> F-7T+) = Use*? + v^Aie"4 | (23a) 



and 



A(D° - ^vr°) = -L {V2A f e W § - Aie^} (236) 



V3 

Having thus determined Ai and A3, we reinstated the phases with i 5x — 63 = 

22 22 

(86 ± 8)° and calculate the branching ratio. Fitting, the branching ratio data deter- 
mined a\ e ^ and a 2 e -^ to lie in the following ranges 13, 

1.13 < ai e// < 1.17, - 0.46 < a 2 e// < -0.42 (24) 

We repeated this procedure for D — > K*tt and Kp decays where the difference of 
isospin phases is known land obtained!, 

K*tt : 1.74 < ai eff < 1.96, - 0.53 < a 2 eff < -0.43 (25a) 

and 

Kp : 1.17 < ai e// < 1.32, - 1.00 < a 2 e// < -0.75 (256) 

In the above calculation we used the measured form factors as much as possible; 
F DK {0),A 1 DK *{0),A 2 DK *{0) and V DK *{0) from Ref. [10], F Dn {0) from Ref. [8] and 
the as-yet-unmeasured A Dp (Q) from the model of Ref. [1]. From Eq. (24) we note 
that the effective a\ and a 2 are very close to the values advocated in Ref. [1] which 
gave rise to the belief that Nq — > 00 limit was relevant to D decays. However the 
value of and a 2 e ^ determined from K*7r and Kp decays belie that belief. 

4. Nature of a^^ and a 2 e ^ 



Up to this stage our discussion has assumed that the effective and a 2 e -^ are 
real . Here we discuss how final-state interactions can be introduced and how they 
render a\ e ^ and a 2 e -^ complex . 
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Figure 3: Inelastic final-state interaction Feynman diagrams 

To illustrate the ideas, we look at the decay D$ + — > 07r + .It involves a single 
isospin amplitude, I = 1, in the final state. Before one introduces inelastic final state 
interactions, the decay amplitude for D s + — > 07r + is given by 

G p 



A{D S + -> 07T+) = -|v cs K>i e// (2m )/.6 • PDs A D ^(m2 ) (26) 



where 



eff _ n J i ^ A) Wn/ L _ .//». X 



ann 



^ " 01 i 1 + V*V*) + a, A^ml) + /J V s Vj) J (27) 

where apart from the nonfactorized contributions we have also absorbed a real anni- 
hilation term in the defenition of an . 

We next unitarize the amplitude through final-state interactions . This renders 
the amplitude complex and one can define a complex, unitarized as follows ( the 
superscript U stands for 'unitarized '): 

A u {Ds+ _» 07r+) = ^£ VcsV:da u,e ff{2m ^ Ue . p^Ds*^ (28) 

For the sake of illustration, we consider coupling of the 07r + channel to the G-even 
K*K channel given by the symmetric state 

| K*K> S = -j= {| K* + K°>+ | K + K°*>} (29) 
In Figure 3 we have shown schematically how the channels are coupled. 



We adopt the K-matrix formalism to unitarize the decay amplitude. The uni- 
tarized decay amplitude, are given in terms of the following coupled equations, 

AU,iK-K )8 ) = {l-&K) { AiK * K)s ) (30) 

where K is a real, symmetric 2x2 matrix with dimensions (GeV)~ 3 ,k. 3 is a diag- 
onal matrix with entries kf and k\ as appropriate P-wave threshold factors for the 
two channels and and A^ K * K 's are the un-unitarized amplitudes. The latter is 
constructed out of the following two amplitudes, 

A(D S + -> K»K* + ) = ^V cs V: d a 2 e "f K (2m K ,)e-p Ds A D ^ K \m 2 K ) (31) 

with 

.//_ \^ B ° 1)nf ■ g i B ^ )nf , a ^f°s B Q ann } 

tt2 ~ a2 { + V s ** {ml) + a 2 V s ** (m%) + a 2 f K V s ** «) ) ^ 

and 

A(D S + - = ^KsK>2 e// /i.* (2m K .)e ■ p Ds F 1 D " K \m 2 K *) (33) 

with 

r {1)n/ r» r {8)n/ n f r ann 

S (1)n/ , 5 (8)n/ and B ann are the analogue of A {1)nf , A (8)nf and A ann of Eq. (27). 
The hatted quantities refer to the channel K° K + . Using the above equations the 
amplitude A^ K * K ' S is written as, 

Al*-«)s = ^V cs V ud *^l {a 2 ^f K A D s K *(ml) + aV ff f K ^ K (m 2 K .)} (35) 

The i^-matrix being real and symmetric is parametrized as follows, 

*=(l "') (36) 



b c 

On carrying out the unitarization through Eq.(30) one obtains 



with 



A U (D S + - 07T+) = ^VcJCfli^'fe ■ PDs (2m 4> )A D ^(ml) (37) 



A I \f2ms 1 



where 



F 



a 2 e// fK Ao D ° K \m? K ) d 2 eff fK. F 1 D s K (m 2 K .) 



and A = det (1 -ik 3 K). 

If the final-state interactions were elastic, b = c = and A = 1 — iki 3 a, we would 
have obtained 



In summary, one can always define an effective, and complex a\ through the 
expression in Eq.(37). What we have shown is how effects such as nonfactorized 
contributions, annihilation and final-state interactions are built into it. A corollary 
of our point of view is that claims of test of factorization by comparing two-body 
hadronic rates to semileptonic rates are really nothing more than determinations of 
| ai u ^f \. 

I wish to acknowledge collaboration with A. B. Santra and F. Ghoddoussi in the 
research described here and thank the Natural Sciences and Engineering Research 
Council of Canada for an Individual Research Grant. 



VI + h 6 a 2 

with 5 = tan~ 1 (aki 3 ), the elastic P-wave 07r scattering phase. 




(39) 
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